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The Holstein-Primakoff representation for spin systems is used to derive expressions with solutions 
that are conjectured to be the zeros of Hermite polynomials H„{x) as n ^ oo. This establishes a 
correspondence between the zeros of the Hermite polynomials and the boundaries of the position 
basis of finite-dimensional Hilbert spaces. 


The Hermite polynomials are prevalent in many fields. 
They can be defined as 

V ) ■ <'> 

In the physics community, they are perhaps best recog¬ 
nized as the Gaussian-weighted eigenfunctions (in posi¬ 
tion representation) of the quantum harmonic oscillator 
(with h = 'm = uj = \,B. convention that will be used for 
the rest of the paper): 
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Associating a spin with a boson c^, Holstein and Pri- 
makoff showed that to satisfy these commutation rela¬ 
tions, the operators can be expressed asQ 


^^ = ctc-^, (4) 

S~^ = \/2S — cfc, and S~ = \/2S — &cc. (5) 


As such, they are orthogonal over the Gaussian- 
weighted whole domain, Hm{x)Hn{x)e~^ dx = 

^/7r2"n!Snm- This last property allows their use in Gaus¬ 
sian quadrature, a useful and popular numerical inte¬ 
gration technique where f{x)dx is approximated as 

fi^j) where Xj are the zeros of Hn{x) and 
f{x) is a well-behaved function. For this and many other 
reasons, an analytic formula for the asymptotic zeros of 
Hermite and other orthogonal polynomials has been a 
subject of much interest [J-Q, especially in the applied 
mathematics community and the field of approximation 
theory. 

In this paper, I examine the position state representa¬ 
tion of the eigenstates of finite dimensional S'-spin sys¬ 
tems, as expressed in the Holstein-Primakoff transforma¬ 
tion. As S' 00 , the system becomes the infinite dimen¬ 
sional harmonic oscillator. This association allows me to 
derive the simple main results presented in eqs [6] and 0 
with solutions that I conjecture become the asymptotic 
zeros of the Hermite polynomials (as n —>■ 00 ). Further¬ 
more, I numerically show that this convergence is rather 
quick and so the expressions can frequently be used, in 
many instances of finite-precision application, as the ef¬ 
fective zeros of Hn{x) with finite n, such as in applica¬ 
tions of Gaussian quadrature. In a more aesthetic sense, 
these results also establish a beautiful correspondence be¬ 
tween the boundaries of equal area partitions of circles 
with radii that are increasing in a certain manner and 
the Hermite polynomial zeros. 

Spin systems are defined by the fundamental commu¬ 
tation relations between operators S^, S'^ and S~ : 


This is a very useful association and has found many 
applications in the condensed matter field’s study of 
many-body spin systems. Each boson excitation repre¬ 
sents the “ladder up” finitesimal excitation away from 
the spin’s extremal S state. The Hilbert space is finite¬ 
dimensional and possesses 25” -|- 1 states {—S,—S + 
I,..., S'}. In fact, considering eq. [5] it is clear that the 
Hilbert space outside this defined space is not even Her- 
mitian. 

Transforming from the Holstein-Primakoff bosonic rep¬ 
resentation to position (and its conjugate momentum) 
space (using the relations = ■^ (<? ~ *p) and c = 

id + ip)) reveals that the trivial Hamiltonian is the 
harmonic oscillator: H = Sz = \{q^ + p^) — (S + ^). 
Moreover, transformation of the S'*" and S~ in eq. [5] 
reveals that the Hilbert space spans the domain = 
p^ + < V4S -I- 1. Just as in the Sz representation, 

25' states all with the same area must exist within this 
domain. Fig. [T] sketches out what they look like for the 
{5 = ^, 5 = 1,5 = |}-spin systems. 

For a particular 5-spin system, the lowest eigenstate 
must have the same sign at all g-basis elements since 
it must be nodeless. On the other hand, the highest 
eigenstate must have n — 1 nodes and so the g-basis ele¬ 
ments must alternate in sign such that the eigenfunction 
passes through zero between them. This latter behav¬ 
ior is sketched in fig. [T] in red by the Hermite polynomial 
Hn{x) denoting the value of the overlying g-basis element 
for the highest eigenstate. 

For 5 —>■ 00 , the Hilbert space becomes infinite¬ 
dimensional and the Hamiltonian becomes that of the 
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FIG. 1. The g-basis representation of a) S = b) S = 1 
and c) 5 = I systems is shown. The radius of the Hilbert 
space’s domain is equal to \/AS + 1 and so grows along with 
the nnmber of allowed basis elements. 


harmonic oscillator defined over (p, q) G with the as- 
sociated eigenfunctions proportional to e~~Hnix). It 
therefore follows that as ^ oo, the boundaries of the 
g-basis elements become the zeros of the Hermite poly¬ 
nomial Hn{x) where n = 2S since the highest eigenstate 
must still have alternating sign with each q-basis element. 

Hermite polynomial zeros xj are real and symmetric 
around a; = 0. To determine these boundary points, the 
25'-dimensional Hilbert space’s circular shape in position 
space can be exploited. For even 25', the area of the all 
the q-basis elements up until the jth boundary (measur¬ 
ing from the origin) is For odd 25, the area is 

7rr^^^. This is illustrated in fig. [H 



FIG. 2. The area of the central a) 2j — 1 or b) 2j q-basis 
elements that approximately determine the jth zero of the 
Hermite polynomial Hnix) for n even and odd respectively 
is shaded in blue. The approximate jth zero is at the right 
boundary of these regions. 


Using simple relations for the area of circle sectors and 
rectangles, it is possible to relate these q-basis element 
areas to Xj] The equation involving the approximate ze¬ 
ros of Hermite polynomials with n even is: 


(2j - l)7r 
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while for odd n it is: 


2j7r 
n -I- 1 


= sin0 -b 6, 


(7) 


where 9 = 2 sin“^ ^ and r = \/2n -I- 1. 

Solving these equations for Xj yields the approximate 
jth zero for the nth Hermite polynomial. The results for 


the zeros of the first 50 Hermite polynomials are com¬ 
pared to the exact zeros in fig. |3l In both cases, eqs. [6] 
and[7]converge to the zeros of the Hermite functions quite 
quickly 0. 



FIG. 3. Exact jth zeros of the Hermite polynomials Hnix) 
for n a) even and b) odd compared to those obtained from 
solving eqs. |5] and [T] 


The finding that the boundaries of equal area par¬ 
titions of growing circles correspond to the asymptotic 
zeros of the Hermite functions appears to be a novel 
one from a search of the literature. It is all the more 
surprising that the origin of this one-to-one correspon- 
dance stems from the Holstein-Primakoff representations 
for finite-dimensional spin systems. Furthermore, on a 
practical level, the apparently rapid convergence of these 






































3 


solutions suggests that they may be useful for more ef¬ 
ficient determination of Hermite polynomial zeros for 
large-dimensional implementations of Gaussian quadra¬ 
ture. 
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